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ABSTRACT 

We present a method for parametric modelling of the physical components of the 
Galaxy's magnetised interstellar medium, simulating the observables, and mapping out 
the likelihood space using a Markov Chain Monte-Carlo analysis. We then demonstrate 
it using total and polarised synchrotron emission data as well as rotation measures of 
extragalactic sources. With these three datasets, we define and study three components 
of the magnetic field: the large-scale coherent field, the small-scale isotropic random 
field, and the ordered field. In this first paper, we use only data along the Galactic 
plane and test a simple 2D logarithmic spiral model for the magnetic field that includes 
a compression and a shearing of the random component giving rise to an ordered 
component. We demonstrate with simulations that the method can indeed constrain 
multiple parameters yielding measures of, for example, the ratios of the magnetic 
field components. Though subject to uncertainties in thermal and cosmic ray electron 
densities and depending on our particular model parametrisation, our preliminary 
analysis shows that the coherent component is a small fraction of the total magnetic 
field and that an ordered component comparable in strength to the isotropic random 
component is required to explain the polarisation fraction of synchrotron emission. We 
outline further work to extend this type of analysis to study the magnetic spiral arm 
structure, the details of the turbulence as well as the 3D structure of the magnetic 
field. 

Key words: ISM: magnetic fields - Galaxy: structure - polarisation - radiation 
mechanisms: general - radio continuum: ISM 



1 INTRODUCTION 

Observations of external galaxies show that the relationship 
between the magnetic field and the ionised gas is far from 
simp le and varies from galaxy to galaxy. (See, e.g., iBeckl 
2009 for a review.) Our own galaxy is more difficult to study 
since we must observe it from within and looking through 
the plane. In this work, we present a method for comparing 
magnetic field models of many parameters to a variety of 
observables. With the tools we describe, we can use more 
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realistic models for the small-scale components and there- 
fore have less need for the simplifying assumptions (such as 
isotropy) that are often made with unknown effect on the 
results. For the first time, we now have coverage of signif- 
icant portions of the Galactic plane in total and polarised 
synchrotron intensity as well as in rotation measures (RM), 
and we show how these complementary datasets are vital 
for disentangling the different components of the Galactic 
magnetic field. 

It is useful to think of the Galactic magnetic field as 
separable into three components referred to as coherent, or- 
dered, and random/tangled/turbulent. These are illustrated 
by the cartoon in Fig.[T] The term "coherent" refers to, e.g., 
a large-scale spiral structure, while the "random" compo- 
nent usually refers to the small-scale component varying in 
three dimensions in both strength and direction. The term 
"ordered" refers to a field where variations simply imply 
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Figure 1. Cartoon illustrating the three components of the mag- 
netic field and how they relate to the three observables of total 
synchrotron intensity (I), polarised synchrotron intensity (PI), 
and Faraday rotation measure (RM). The variance of these is 
also of interest, as shown in the case of ctrm- (Note that 
the situation is the same for dust emission, which is also polarised 
perpendicular to the magnetic field due to grain alignment.) 



sign reversals; the ordered component points along a com- 
mon, ordered axis but simply changes direction on small 
scales, perhaps stochastically. (There are different uses of 
the word "ordered" in the literature. It can be thought of as 
an anisotropic random component, but in this paper, we will 
distinguish between an isotropic random component and an 
ordered component. "Ordered" is also sometimes used to 
refer to the combination of what we in this paper call the 
ordered plus coherent field, but we find that usage confusing 
and prefer to think of them as three distinct components.) 

The coherent fields are assumed to probe processes 
such as galactic-scale dynamos, while the small-scale ran- 
dom component probes turbulent processes in the ISM. The 
ordered component is thought to result from the effects of 
larger scale shearing dynamics, compression, etc. on the ran- 
dom component. This comp onent is often neglected for sim- 
plicity (see, e.g.. lBurr]|l966l ). 

In external galaxies, the synchrotron emission is com- 
pared to the thermal emission to study how the magnetic 
fields correlate with the ionised gas in the galaxy. The arms 
seen in total intensity tend to follow the arms seen in thermal 
emission, but the same is not true of the polarised intensity. 
This has two implications. Firstly, the random component is 
often significant, perhaps dominant, in the spiral arms, thus 



lowering the polarisation fraction. Secondly, the ordered and 
coherent fields may be strongest in "magnetic arms" distinct 
from the spiral arms traced by the diffuse ionised gas (DIG) 
component. As reviewed by iBeckl (|2009l ). NGC 6946 is an 
example where the magnetic arms appear to be between the 
arms traced by thermal emission, while M 51 appears to 
have magnetic arms on the edge of the spiral arms. 

Though we cannot assume anything about how the 
magnetic field components correlate with the DIG, Fig. [2] 
shows that the synchrotron emission has clear step features 
that have b ee n thought to correspond to arm tangents (e.g., 
Mills 1959] or iBeuermann et all TT985). The rotation mea- 
sures also indicate field reversals that may be related. Learn- 
ing how the field varies in strength, direction, and coherence 
across a spiral arm traced by diffuse gas will inform theories 
of how the fields are generated and maintained in the dy- 
namic environment of the magnetised interstellar medium 
(MIM). 

In this work, we take a first step toward disentangling 
these components in our own Galaxy. Since we cannot look 
down upon it from above to determine where the magnetic 
arms lie relative to the DIG spiral arms, we are limited to 
looking through the plane. Our ongoing project is to see 
whether the profiles in total and polarised synchrotron emis- 
sion can distinguish arm ridges in the different components, 
but for this work, we assume that they all peak in the same 
ridges. We do not, however, constrain those ridges to lie 
where the DIG spiral arm ridges are thought to be, but in- 
stead use the data to constrain the orientation of the mag- 
netic spiral arms. We leave it to later work to determine if 
the data can distinguish models where in addition to inde- 
pendent magnetic arms, the peaks in the three components 
do not coincide. 

But there are several difficulties in interpreting the 
emission profile of synchrotron total intensity along the 
Galactic plane. Firstly, the random component introduces 
a sort of "galactic variance" (analogous in a limited sense to 
the cosmic variance in cosmic microwave background stud- 
ies). Nearby random features will perturb the profile mak- 
ing it difficult to separate the field components accurately. 
Secondly, the plane contains a strong emission component 
from thermal bremsstrahlung, commonly known as free-free 
emission. This component is difficult to separate from the 
non-thermal emission we wish to study and has a different 
spatial distribution. 

Using the polarised synchrotron emission in addition 
helps to resolve some of these problems, as the isotropic 
random component does not on average contribute to the 
polarised intensity. The addition of rotation measures fur- 
ther helps to break the degeneracies in the parameter space, 
as they are only dependent on the coherent component. The 
three field components can then be studied by comparing 
these three observables. 

It is only recently that we have enough data to study 
the three components of the magnetic field and perhaps to 
separate them. We can look at the continuum emission and 
compare the total intensity to polarised intensity as well to 
the Faraday rotation measures. These three contain infor- 
mation about the relative strengths of the coherent, ordered, 
and random components. We now have synchrotron emission 
maps over the full sky at several bands including polarisa- 
tion, and we have unprecedented coverage of a large part 
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Figure 2. The three primary tools we use for studying the mag- 
netic field in the plane of the ga l axy a re the synchrotron total in- 
tensity (top) from lHaslam et al.l( 1982), the synchrotron polarised 
intensity (middle) from Hinshaw (2009), and the RM data (bot- 
tom) from lBrown et all j2003l , 120071) averaged into roughly 6° 
bins. The field structures that we can perhaps infer from the pro- 
file of these data along the plane are indicated. The vertical lines 
mark interesting sight-lines (solid for positive RM, dashed 
for negative), also shown on Fig. [4] showing what may be tan- 
gents to arm features. Also shown in orange is the estimate for 
the free-free contamination which has been subtracted from the 
total emission at 408 MHz (dotted green line); see § 12,11 



of the Galactic plane in RMs of extragalactic sources. The 
combination of these datasets allows us to probe the three- 
dimensional distribution of magnetic fields and ionised gas 
in the Galaxy. 

There has been much previous work to mode l the 
Galactic magnet ic field , fro m iBeuermann et al.l ll 19851) and 
iBroadbent et alJ ill990l) to lHan et alj ll2006l), iPage et al.l 
ll2007l). IMiville-Deschenes et all (|200S|), Isun et alj i|2008l ). 
iJansson et all (I2009T ). and lOrlando et al.l (I2009T ). Most of 
these works use only one or two datasets and are limited 
by simplifying a ssumptions about the n ature of the irregu- 
lar ^cojruoonmt^^roadben^eF^O (1990) (with more details 
in iBroadbentl Il989l ) consider a detailed galaxy model and 
have remarkable success in reproducing the profile of syn- 
chrotron emission along the plane. That work is unusual 
in including an analytic method to model the anisotropy 
in the small-scale random component of the MIM. While 
it seems unlikely that the data constrain all of the dozens 
of parameters in their model (particularly considering the 



above-mentioned problems of thermal emission separation 
and galactic variance), we consider it an interesting place 
to start and ad opt severa l aspec ts of their model. The more 
recent work of ISun et alJ (2008) uses the three observables 
together to present a 3D model of the galaxy, but it does 
not explore the parameter space systematically and consid- 
ers only what we define as the coherent and isotropic ran- 
dom m agnetic field components. The work of IJansson et al.l 
(2009) does a systematic search using MCMC, but they do 
not explicitly treat the random (or ordered) components and 
use only two of the three datasets we consider vital. 

Our principle aim in this work is to describe a method 
to model the physical components of the MIM, to simulate 
the observables, and to map out the likelihood space using a 
Markov Chain Monte-Carlo (MCMC) analysis. The problem 
is a large and complicated one, and in this first paper, we 
focus on a 2-dimensional analysis on the Galactic plane only 
(seen as a ID profile as in Fig. [2} , though the methods can 
certainly be extended to model the Galaxy in 3D. We then 
demonstrate the method's utility by using the available data 
to constrain the relative strengths of the three components of 
the magnetic field, the coherent, random, and ordered fields, 
and discuss how to extend the analysis using additional data. 



2 OBSERVATIONS 

2.1 Synchrotron emission in total intensity 

The 408 MHz full-sky map bv lHaslam et alj l|l982l) primar- 
ily consists of synchrotron emission. (This map is available 
in HEALPbQ format from the LAMBDA websitfl) The in- 
tensity gives information about the magnetic field (B) com- 
ponent perpendicular to the line of sight: 



z (v,p) oc 



LOS 



JcR E (x)B ± (x) P ^ 1 d/ 



(1) 



where / is the specific intensity and Jcre(x) is the density 
of cosmic ray electrons (CREs) described explicitly in § 13.51 
The total intensity, or Stokes I, in a given observing beam 
is then 1= j IdQ. 

The index on B depends on the spectral index, p, of the 
cosmic ray electron distribution, whic h we assume to be a 
power law. We follow the notation of iRvbicki fc Lightmanl 
(1 19791 ) and define the number density of particles in the 
range 7 to 7 + d"/ (where 7 is the Lorentz factor such 
that E = jm e c 2 ) as ^(7)^7 oc r y~ p d r y. Observations show 
that p w 3 between our two synchrotron o bserving frequen- 
cie s of 408 MHz and 23 GHz ( see, e.g., iFinkbeinerl (20041 
or IMiville-D eschene s et al]|2008l ). in which case dl oc B 2 . 
It also implies a power law in synchrotron flux density of 
S oc j/~( p-1 )/ 2 = v~ x or brightness temperature of T oc v~ A . 
In reality, the spectral energy distribution of electrons will 
vary with position in the galaxy, but we will approximate it 
as a power law with a constant p = 3. The impact of this 
choice is discussed further in § 16.61 

The synchrotron total intensity is sensitive to all com- 
ponents of the magnetic field. The coherent component, ran- 
dom component, and ordered component all contribute ac- 
cording to their projection perpendicular to the line of sight. 



1 http://healpix.jpl.nasa.gov/ 

2 http://lambda.gsfc.nasa.gov/ 
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The profile of the synchrotron emission along the Galactic 
plane, shown in Fig. [2] has distinct steps, exactly as if the 
emission increases as the observer looks along a spiral arm. 
Unless this is coincidental, it therefore seems that either the 
magnetic field structure or the CRE density must follow to 
some degree the DIG spiral structure. 

The 408 MHz map includes thermal bremsstrahlung 
(free-free) emission in a very narrow region about the Galac- 
tic plan e. (Other available radio surveys include the 2. 3 GHz 
data of I Jonas et al.lll998 | or the 1.4 GHz survey of iReichl 
1982 and lReich fc Reich| [l986 Higher frequencies, however, 
will have more free-free contamination due to its flatter 
spectral index.) Since it is along the plane that we are 
interested in the synchrotron emission, we must subtract 
the thermal emission. We use the Wilkinson Microwave 
Anisotropy Probe ( WMAP) foreground maps generated 
from the maximum entropy method (MEM) described in 
iHinshaw etUI lj2007l ). The WMAP MEM free-free map for 
the Ka-band (33 GHz) contains the best estimate of the 
free-free distribution that we have for the Galactic plane 
region. f (|de Oliveira-Costa et all [2008). for example, pro- 
duce a synchrotron map using a combination of 11 different 
frequencies, but their principal components analysis results 
in the 408 MHz profile along the plane being virtually un- 
changed, indicating that the free-free component has not 
been subtracted.) The MEM maps are also available from 
the LAMBDA website. 

The data are generated by starting with the full sky 
HEALPix maps at JV a id e = 512 (pixel size of 6.9 arcmin) for 
the Haslam and the WMAP Ka-band MEM free-free compo- 
nent. Each map is then smoothed to an effective beamwidth 
of 3° FWHM and the maps downgraded to N sidc = 128 
(pixel size of 27.5 arcmin). From the Haslam data, we then 
subtract the free-free component extrapolated from the Ka- 
band assuming a pow er law dependence of v~ 21 (see, e.g., 
iDickinson et al.ll2003h . The 512-pixel (4Af sidc ) slice along the 
plane is then extracted, further smoothed in 1-dimension 
by boxcar-averaging by 8 pixels, and the result is then 
downgraded to the 64 bins whose centres coincide with the 
-Wside = 16 pixels along the plane. The result is a profile with 
a resolution of ~ 6° in longitude and three in latitude. 

The reason for this processing is to minimize the effects 
of the inaccurate thermal emission separation. The WMAP 
MEM free-free component is generated by an analysis in- 
cluding only three foregrounds, namely synchrotron, free- 
free, and thermal dust. The results are therefore contami- 
nated by the anomalous dust component. There is no other 
reliable indication of the thermal emission on the plane, how- 
ever, (Ha, for example, is absorbed by dust on the plane and 
is only a reliable tracer of free- free at high latitudes) and the 
free-free emission is strong enough on the plane relative to 
the other components that the separation, while imperfect, 
is probably suffic ient. The free- free lat itude profile is very 
narrow (see e.g.. IDickinson et all 120031 ) compared to that 
of the synchrotron emission. Smoothing the maps first then 
leaves the profile of synchrotron emission along the plane 
essentially unchanged but reduces the relative amount of 
free- free emission. 

In future work, we will investigate other ways to sep- 
arate the thermal emission on the pla ne, for example b y 
using the newer results from WMAP in iGold et all ll2009h . 
which include a pixel-based MCMC separation method. In 



that work, the authors estimate that the spinning dust frac- 
tion is at most ~ 15 — 20 per cent of the emission in the 
Ka-band. It is hard to predict the fraction that ends up in 
the free-free template (in either the MCMC or MEM analy- 
ses) when the spinning dust component is ignored, since this 
depends on whether its spectral behaviour better matches 
the free-free or the synchrotron. Note that the remaining 
uncertainties due to the inaccurate thermal emission sepa- 
ration and subtraction are unlikely to be as significant as 
the galactic variance described in § [1] itself. 

Because the Galactic centre region is likely complicated 
and currently mysterious, we exclude the four nearest pixels 
corresponding to roughly 10° either side of £ = 0. This makes 
the analysis insensitive to what happens in the innermost 
~2 kpc of the galaxy and allows us to focus on the step- 
features as well as the general profile. 

2.2 Polarised intensity 

Polarised synchrotron emission depends on the cosmic ray 
power law spectral index, and is at most a fraction n of the 
total intensity, where 



n 



p+1 
p + 7/3 



= 0.75 



(2) 



(for p = 3) in the case tha t the magnetic field is uniform 
ijRvbicki fc Lightmanlfl979T ). The observed degree of polari- 
sation compared to the total intensity therefore gives infor- 
mation about how ordered the magnetic field is perpendicu- 
lar to the line of sight. Averaging over an isotropic random 
component produces, on average, no polarised intensity, as 
emission polarised in perpendicular directions cancels. An 
ordered or anisotropic random component, however, will still 
add to the polarised intensity, since the polarisation angle 
depends only on the orientation rather than the direction of 

the field. 

We use the IHinshaw! (|2009l ) five-year WMAP 23 GHz 
K-band map of polarised intensity and smooth it to extract 
the plane as described above for the Haslam data. Likewise, 
we ignore the four pixels toward the Galactic centre. This 
dataset is currently the best polarised synchrotron map in 
the high frequency regime where Faraday effects are negli- 
gible. 

2.3 Faraday Rotation measure 

The polarisation angle of an electromagnetic wave rotates 
when propagating through a magnetised plasma. The rota- 
tion angle changes as RM\ 2 , where the Faraday rotation 
measure (RM) is defined as the the line-of-sight (LOS) in- 
tegral from a source at a distance D: 



RM oc / n e B\\Al 
J D 



(3) 



where n e is the number density of thermal electrons. A pos- 
itive RM means that Bm points toward the observer. 

The RM can be measured for both pulsars within the 
galaxy as well as external sources. These data complement 
the synchrotron emission, since the RM is sensitive to direc- 
tion as well as orientation and to the parallel component 
rather than the perpendicular. Therefore, averaged RMs 
trace only the coherent component of the magnetic field. If 
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the thermal electron density is uncorrelated with the mag- 
netic field, then any random component or even ordered 
component will have a null effect on average, though it will 
naturally give rise to a variance in the manner of a random 
walk. This means that these data can help break degenera- 
cies that arise in modelling the synchrotron emission due to 
the coherent versus random field parameters. 

We are therefore probing three components of the mag- 
netic field using three main datasets. There is a further 
source of uncertainty, however, namely the distribution of 
thermal electrons in the galaxy, discussed in § 13.41 As men- 
tioned above, the relationship between the DIG and the 
magnetic field structure is unclear. 

The data we use are the RMs for the extragalactic 
sources in the Can adian Galactic Plane Survey (CGPS) of 
iBrown et all d2003Tl and in the Southern Galactic Plane Sur- 
vey (SGPS) of lBrown et"afl (|2007T l. In the CGPS, there are 
a total of 380 sources within 5° of the plane in the range 
82° sC £ sC 146°. In the SGPS, there are a total of 148 sources 
within 1.5° of the plane in the range 253° ^ i ^ 356°. Fol- 
lowing Brown et al., we exclude from our analysis the 30 
sources in the region 270° < t < 280°, where there is evi- 
dence for anomalously low RMs due to a local feature. To 
maintain consistency with the other datasets, we only use 
the 151 CGPS sources within 1.5° of the plane. That leaves 
a total of 269 sources used in our analysis covering a bit less 
than half of the Galactic plane. The sources are then aver- 
aged into the bins corresponding to those used for the other 
datasets. 

There are additional RM data for galactic pul- 
sars that may in future be added to our analysis. 
Because each pulsar is at a different LOS distance 
through the galaxy model, to include them would 
complicate the analysis and add considerable pro- 
cessing time. For this work, then, we use only ex- 
tragalactic sources that give a full LOS through the 
galaxy. 



3 GALAXY MODELS 

We must model each of the relevant components of the MIM, 
namely the magnetic field (coherent, random, and ordered), 
the thermal electron spatial distribution, and the cosmic ray 
electron spatial and spectral distribution. 



3.1 Coherent magnetic field 

We model the coherent magnetic field beginning with a sim- 
ple axisymmetric spiral that defines the direction of the field. 
In Sun-centric coordinates, then, the field direction is 

B = sin(6» p +(/>-£)&- cos(6> p + <j> - l)y 

where cj> is the azimuthal angle in the polar coordinate sys- 
tem with the Sun at the origin, while i is the Galactic lon- 
gitude. We are working only in the plane, and the vertical 
component is always zero. (We use this form since we will 
integrate along the LOS from the Sun as the origin.) 

This defines a spiral field with a pitch angle of 9 P = 
— 11.5. This angle is roughly that used in the thermal 
electron density model described in § 13.41 Estimates of 
9 P in the literature vary, however. lHan 8z Qiaol 



(jl994h get dp = -8°.2±0°.5 b ased on RM s urveys, 
in rough agreement with the iHeilesI (| 19961 ) result 
based on starlight polarisatio n, 8 P = —7°. 2 ± 4°.l. 
IMiville-Deschenes et al.l (| 20081 1 find a value of 8.5° 
based on WMAP polar ised synchrotron emission , 
while IPage et ail (|2007h and Ijansson et al.l (2009) 
find best-fit pitch angles of as much as 35°. Jansson 
et al., however point out that this is highly model- 
dependent. In later work, we will allow this param- 
eter to vary to see if it can be constrained better, 
and to investigate whether the magnetic spiral arms 
follo w the matter s piral arms. 

iBroadbent et all (|l990h defined the coherent field am- 
plitude as: 

B(r) = B (l-exp(-r 2 /i?2))(exp(-r- 2 /i?o)+exp(-r 4 /i?^)) 

(4) 

where r is the Galacto-centric radius, and the three scale 
radii define the shape. Outside of the gala ctic centre r egion , 
this is a simple expo nential drop as used in Han et all (I2006T ) 
to fit RM data or in IMiville-Deschenes et al.l (|2008l ~for po- 
larised emission. Because of its success predicting the syn- 
chrotron total in t ensity profile, even in the inner galaxy, in 
IBroadbent et alJ (119901'), we start with this mo del. But we 
note firstly that iMiville- Deschenes et~ai1 (|2008h find a con- 
stant amplitude fits the WMAP data better, and secondly, 
this profile is degenerate with the cosmic ray density profile 
described below in § 13.51 (In fact, Broadbent et al. use a con- 
stant CRE density, so their profile was essentially accounting 
for any radial variations in both components.) 

Constraining this profile is beyond the scope of this pa- 
per, as it will require an additional observable such as dust 
emission to break the degeneracy. The default values we use 
are given in Tabled] though as described in S I2.1I we exclude 
the innermost longitudes from the analysis, which means the 
analysis is not sensitive to Ri and R2 ■ The value for Ro is de- 
generate with the h r parameter controlling the cosmic ray 
profile defined in § 13.51 also an exponential disc. The de- 
fault values in Table [1] were set by eye to approximate the 
synchrotron profile so that we could explore the parameters 
determining the random versus ordered components. In a 
later paper, we will explore the radial profile parameters in 
more detail. 

This azimuthally symmetric field is then amplified along 
spiral arm ridges (defined parallel to the field direction) as 
described in § 13.31 Note that those ridges are not constrained 
to lie along the DIG spiral arm ridges, and each magnetic 
arm has a corresponding amplitude parameter, a n , so that 
we can fit for the strength of each arm independently. The 
parameters are summarised in Table Q] and the radial profile 
shown in Fig. [3] 



3.2 Random magnetic field 

The polarisation of synchrotron emission in the spiral arms 

of external gal axies is generally low (a few percent according 

to lBeckll2009l ). implying that a large fraction of the total 

magnetic field is due to an isotropic random component. 

The random component we simulate begins with a 

Gaussian Random Field (GRF) simulated in Fourier space 

(k is the wavenumber, 1/A) with a variance given by 
| fe |(a-2)/2 ) where 

a is the spectral index of the power law de- 
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Figure 3. An example of the radial profile of the coherent field 
and of the compression. (Note that the values used are just for il- 
lustration and not those found to fit the data.) The uncompressed 
profile is shown in red, while the dotted black line includes the 
compression factor, p c (d), showing example spiral arms. The bot- 
tom axis is the radius from the Galactic centre, while the top axis 
is the corresponding longitude if viewed from the perpendicular 
direction at a distance of = 8.5 kpc. 



scribing the t otal magnetic ene rgy in one dimension, P(k) oc 
k a , following lHan et al] i|2004h . For a Kolmogorov-type sim- 
ulation that perhaps matches the small- s cale (1/fc < lpc) 
fluctuations, a = —5/3, while lHan et al] ()2004r ) find an in- 
dex of a = —0.37 at larger (1/fc Z lkpc) scales. 

This field in fc-space is then Fourier transformed to real 
space and amplified along the ridges of magnetic spiral arms 
by a parametrized compression described in the next section. 

In this paper, we are simulating the entire galaxy, 
though only in two dimensions. In § [5] we describe how the 
resolution of the random component affects the results. In 
short, we find that the analysis of the whole plane at low res- 
olution is not strongly sensitive to many parameters of the 
GRF. We use the spectral index for larger scales, a = —0.37, 
a box size 40 kpc long on each side, and 512 bins for a bin 
size of 80 pc. This is very low resolution and not much better 
than a single scale random component. But as described in 
§ we find that increasing the resolution does not change 
the results significantly. 

For later papers, we will analyse small regions of the 
sky at high resolution (e.g., pixels of 0.5 arcmin) and look 
in detail at the properties of the turbulent component, but 
for this paper, those details are not constrained. 



3.3 Compression 

As discussed above, in some external galaxies, the magnetic 
field appears enhanced in the spiral arms as traced by the 
DIG, while in others, there appear to be separate "mag- 
netic arms", sometimes between the DIG arms or showing 
a different (even varying) pitch angle. We therefore define a 
logarithmic spiral arm model for such ridges parallel to the 
coherent field direction and allow the orientation of the spi- 
ral to vary to fit the data. In addition to the spiral arms, we 
also include an annulus at 5 kpc from the Galactic centre, 
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Figure 4. An example of a coherent spiral arm magnetic field 
model (no random component shown) compared to NE2001 spiral 
model for the thermal electron density shown as overlaid contours. 
Interesting sight-lines are over-plotted corresponding to the ver- 
tical lines in Fig. [2] The peaks in the RMs correspond to looking 
roughly tangcntially to a spiral arm. (Note that the contours do 
not make clear the high-density "molecular ring" between roughly 
2 and 5 kpc of the Galactic centre.) For comparison with Table [JJ 
the Perseus arm corresponds to amplitude ciq, Sag-Carina with 
a\, Scut-Crux a.2, Norma (13, and the ring 04. (As discussed in 
§ I3.ll the parameters used in this example are not necessarily the 
defaults in that table.) 



roughly coinciding with the "molecular ring" feature in the 
NE2001 model. 

We define a four-arm model in this paper, and will test 
other models in a later work. Each arm is then described by 
a curve 



r(0) = R 3 exp [(0 O n - <t>)/P] 



(5) 



where (9 = 1/ tan(0 p ), P being the same pitch angle defined 
above for the coherent field direction. The angle 0on is the 
reference angle for the n th spiral arm, always 7r/2 away from 
its neighbours in the case of four arms. 

The arm ridges (and ring) are then defined by a mod- 
ulation that peaks along the ridge and reaches a minimum 
be tween the arms. The co mpression profile we use is based 
on iBroadbent et al' I (|1990T ). who modelled a magnetic field 
component that is enhanced in the arm by a compression 
factor p c , 



p c = C exp(-(d/do) ) + 1 



(6) 



where d is the distance to the nearest of the four arm ridges 
as measured along the line passing through the Galactic cen- 
tre. 

This implies a field that is unaltered between the arms 
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but amplified in the arms according to a Gaussian pro- 
file, when seen though a radial cross-section, with a width 
of do- This amplitude enhancement profile applies to both 
the isotropic random component and the coherent compo- 
nent. This Gaussian modulation may not accurately repre- 
sent what may in reality be a shock front with a different 
profile, but it is meant simply as a first order approximation 
to an enhancement parallel to the shock plane. (Though we 
refer to the "shock plane", note that p c is not the density 
contrast across the shock, as often used, but simply the to- 
tal Gaussian density enhancement as compared to the un- 
shocked inter-arm region.) An example of the compressed 
coherent field viewed from above is shown in Fig. [4] 

To simulate an ordered component, the random field is 
then also stretched along the plane defined by the arm ridge. 
This can be thought of as an anisotropic random component, 
or as the addition of an ordered component (in addition to 
the coherent and isotropic random components) whose ori- 
entation is always along the spiral arm but whose specific 
direction changes stochastically. Rather than using the nu- 
merical approximation given in Broadbent et al., we apply 
this compression explicitly to our simulated random compo- 
nent. Namely, we use 

Bi rrc g = /5 c BgRF + /ord(pc — 1)Bqpp = Bi so + B ord (7) 

where B^j is the component of the random field projected 
onto the shock plane, 

-Rproj _ -p, B coh ■ B G RF , , 

B GRF -B coh ^— p , (8) 

and / or( j sets the ratio of the ordered field to the amplified 
but isotropic random component in the arms. (Note that 
/ or d is not the ratio itself, since that changes from the arm 
to inter-arm regions.) The first term in Eq. [7] is then the 
isotropic (but no longer homogeneous) random component 
while the second term is the ordered component. The shock 
plane is assumed to be perpendicular to the Galactic plane 
and parallel to the spiral arm traced by the coherent com- 
ponent. 

The form of this parametrization means that in the 
inter-arm regions, the ordered component goes to zero, while 
the isotropic random component is not amplified but non- 
zero. The combination of Brms and Co determines the am- 
plification in the arm relative to the inter-arm region. In 
particular: 

/r2 \i/2 _ |(-Bgrf) 1/2 = -Brms, as p c -> 1 (interarm) 

\ {Co + 1)5rms, as p c ->• Co + 1 (ridge). 

(9) 

In other words, in the inter-arm region, d is large, p c goes 
to 1 (see Eq.[6]). On the ridge, d goes to zero, p c reaches its 
maximum at Co + 1, and the random component is amplified 
by Co + 1. In between, the amplification is falling off with 
a Gaussian profile of width do- The profile of the ordered 
component, proportional to p c — 1, is the same Gaussian but 
without the offset, i.e., it goes to zero between the arms: 

(-Bord) ^ = 

!0, as p c — > 1 (interarm) 

/ordCo ((Bgrf) 2 ) 172 as p c -> Co + 1 (ridge) 

(10) 



Note that this model approximates the case where the or- 
dered component is due to the compression wave in the arm, 
but it does not represent the possibility that the ordered 
component might arise simply due to differential rotation. 

The ratio of ordered to random is then zero in the inter- 
arm regions; along the ridge, the ratio is 

W*" , Co [2 

(B? )V 2 ~ /ord C + lV3 (LL> 

\ ISO/ 

because ((B G ^) 2 )= 2/3 <B GRF ). 

When looking down a spiral arm, the ordered compo- 
nent does not contribute to the synchrotron emission. The 
impact of this change on total intensity is then seen in the 
relative amount of emission seen looking down an arm rela- 
tive to elsewhere, i.e., in the shape the emission steps seen 
in the longitude profiles in Fig. [2] Since the isotropic compo- 
nent does not contribute to polarisation, the ordered com- 
ponent, along with the coherent field, also determines the 
polarised emission profile. 

The same effect on the synchrotron can be obtained 
with an isotropic compression of the random component 
combined with a stronger coherent component. Amplifying 
the random component along the arms essentially adds an 
ordered component, which as described above, is indistin- 
guishable from a coherent component in synchrotron. The 
two cases can be distinguished, however, by looking at the 
rotation measures as well, where a too-large coherent com- 
ponent will overpredict RM. 

3.4 Thermal electrons 

For the thermal electrons, we consider three possibilities: 

• a physically unrealistic constant density; 

• the NE2001 model; or 

• a parametrized spiral arm structure analogous to the 
coherent magnetic field component. 

The thermal electro n density model t hat ha s been most 
widely used is that of ICordes fc Laziol (|2002l ). known as 
NE2001. Based on observations of pulsars at known dis- 
tances, this model relates the thermal electron density, n e , 
to the dispersion measure (DM) of a give source and its 
distance, D: 

DM = f n e dl (12) 
Jo 

Though this model was intended for determining pulsar 
distances, it is also used for modelling the magnetic field 
through RM, as it is the best existing predictor of the ther- 
mal electron density. 

The constant model is useful as a comparison to deter- 
mine whether or not we can distinguish either of the others 
from the null hypothesis. We have compared the results of 
using a constant density to those using the NE2001 model 
and find that the resulting data profiles are largely indistin- 
guishable. The models and RM data are shown in Fig. [5l 
with the corresponding parameters that were fitted listed in 
Tabled (The synchrotron profiles, not shown, are effectively 
indistinguishable.) The curves are essentially the same, and 
though the individual arm amplitudes change, the effect on 
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Uniform n =0.1 cn 
NE2001 " 
CGPS & SGPS data 




Ion [cleg] 

Figure 5. A comparison of the rotation measures from two mod- 
els fitted to the CGPS and SGPS data, in the one case with a con- 
stant thermal electron density of n c = 0.1 and in the other case 
with the NE2001 prediction. The best-fit parameters are given in 
Table [2] (The dotted lines indicate the galactic variance due to 
the random magnetic field component.) 



the synchrotron emission would be very difficult to distin- 
guish. 

In this work, we use the NE2001 thermal electron den- 
sity model unless otherwise specified. Though it is not likely 
to be correct in detail, it is at least based on observations 
and provides the best current estimate. 

The possibility of denning our own model introduces 
quite a few more parameters to the problem, and those are 
unlikely to be constrained with the data we are using. In 
the future, however, we can include dispersion measures that 
depend only on the thermal electron density (and are cal- 
culated in the integration code we use) and additional data 
about the distribution from Ha emission. The tools we use 
can be extended to include galactic pulsars, comparing that 
additional RM and DM data to better map out the ther- 
mal electron distribution near the Sun. As each pulsar is at 
a different distance, this complicates the LOS integration 
and the likelihood exploration, but in future work, we can 
explore this idea further. 



3.5 Cosmic-ray electrons 

The spatial and spectral distribution of cosmic rays in the 
Galaxy is thought to be fairly smooth u ntil e nergies above 
100 GeV (see, e.g., IStrong et all |2004| and IStrong et all 
l2007h . But it is not clear how smooth the distribution is, and 
local measurements are obviously affected by local sources, 
so it is difficult to determine the average distribution over 
the galaxy. 

For the purposes of constraining the structure of the 
magnetic field, however, we start simple. W e adopt an 
expon ential disc mode l such as that used by iPage et al.l 
J2007fi and ISun et all d200ct ) (motivated by the work of 
iDrimmel fc Spergelll200li r namely 



■/cre(V, z) = JcRE.e exp(-(r - R®)/hr) 



(13) 



where is the Galactocentric radius of the Solar system. 
(In the current analysis on the plane only, we neglect any 
modulation with height.) The normalisation is uncertain, 



but as a first estimate, we r ead off the value at 10 GeV from 
Fig. 4 of lStrong et all (|2007l ) of roughly E 3 J(E = 10 GeV) « 
250 GeV 2 m~ 2 s _1 sr _1 . We then use the value of JcRE.e = 
0.25 { GeV m 2 s s r) 1 to set the co smic ray normalisation. 
For comparison, ISun et all l|2008h use a value of Jcre,® = 
0.4 (GeVm 2 ssr) _1 . 

For this work, we also adopt a simple power law energy 
distribution with the index p — 3. 

We should emphasize that this model is overly simplis- 
tic in assuming a smooth spatial distribution and a simple 
power law spectrum that is the same throughout the Galaxy. 
Neither of these assumptions is likely to be true in reality. 
What we are constraining is then the spatial distribution of 
the product ncRE^? 2 - The degeneracy between spatial vari- 
ations in the two independently is only broken for the coher- 
ent field using the rotation measure data. In a later work, 
we will see how much the two can be separated by adding 
thermal emission from dust to the analysis. 



4 MODEL SELECTION METHOD 

The problem is then to parametrize the Galactic magnetic 
field (as well as the the thermal and cosmic ray electrons 
if possible) and find the model parameters that best fit the 
data. There are, then, a large number of parameters, and 
a brute-force approach like a simple grid search would re- 
quire prohibitive amounts of computing time. Instead, we 
use the far more efficient and flexible method of a Monte 
Carlo Markov Chain (MCMC) analysis. 



4.1 Simulation with HAMMURABI 

The likelihood function we use computes a simple x 2 from 
a comparison of the model to the data along the Galac- 
tic plane as shown in Fig. [2] For the given values of the 
input parameters, we generate the observed emission and 
RM at each pixel by performing a line-of-sight integration 
through a simu l ated g alaxy using the HAMMURABI code of 
IWaelkens et~aH (|2009l ). 

The HAMMURABI code is designed to compute observ- 
ables such as synchrotron emission and dust emission in full 
Stokes parameters while taking into account Faraday rota- 
tion and depolarisation effects. It performs a line-of-sight in- 
tegration through a galaxy simulation, but uniquely, it can 
refine the integration resolution as the distance increases in 
order to maintain a roughly constant physical cell size. This 
is crucial for simulating effects such as beam depolarisation. 
We do not take full advantage of this code in this current 
work, but we will later use it to investigate more thoroughly 
the properties of the turbulent ISM and its effect on the 
observables. 

For this analysis, we simply use one observing "shell" 



3 This is sometimes expressed as a spatial number density nor- 
malisation. The spatial number density (assuming relativistic 
electrons with v ~ c) is then N(E) = J(E) X 4i X 1/c. Of- 
ten, what is given is a normalisation C(r, z) such that JV(7)d7 = 

C{r,z)^/~ v d') where 7 is the Lorentz factor (E = ■ym e c 2 ). There- 

, 3 



fore, C e = AT( 7 ) 7 p 



icr 



assuming p 



J X 
= 3. 



x m a c^ x 



10 GeV ' 



4 x 
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Parameter 


Default 


Equation 


Description 






Coherent magnetic field: B co h(r, </>) = Bo(r)a 


nPc(d) 


B 
Rs 
Op 

4>o 

Rmol 

^max 
narms 
^max 

Alarms 

a n 

do 

Ro 
Ri 
R 2 
Co 


1 A*G 
7.1 kpc 
-11.5 

5. kpc 
20 kpc 
15 kpc 

4 

0.3 kpc 

20 kpc 
3 kpc 
0.5 kpc 
1 


see p 

/3 = 1/ tan(fjp) and 
r(</>) = R a cxp [(4,o - 
See pitch 

See Co 

B coh {r,<j)) = B (r)a n p c (d(r, <f>)) 
p c (d) = Coexp(-(d/d ) 2 ) + l 

B (r) = B (l - cxp(-r 2 //?|))(cxp(-r 2 /^) + cxp(-r 4 /i?f )) 
see Ro 
see Ro 

See do. Co = if r > R^ 1 ™ 


Global amplitude normalisation 

Scale radius of spiral. 

P is the pitch angle of spiral 

(r(4>) gives arm radius at given azimuth) 

Angle representing rotation of spiral around 

axis through Galactic poles 

Radius of "molecular ring" 

Maximum radius, beyond which \B \ = 

Maximum radius, beyond which Co = (see 

Co). 

Number of spiral arms 

Amplitude modulation and direction for a 
given arm, n. 

Defines width of arm for density enhancement, 
p c ; d is the distance to the nearest arm in kpc 
Outer radial profile parameter. See Fig. [3] 
Inner radial profile parameter 

Peak density contrast. Sometimes Co oc a n \ 
see §H2] 


Random and ordered magnetic field 


a 

D co 

firms 
ford 


-0.37 
1 kpc 


P B (k) = <S ran (fc) 2 ) oc k a 

Bran(fc) = for k < l/D co 
Brms = (B 2 an (x)) 1//2 

B ran — PcB r an + /ord (fic l)Bp r oj = Bi so + B orc i 


Power law spectral index of initial GRF; 
default from Irian et alj l|2004f) in ID. Kol- 
mogorov value —5/3. 

Cutoff maximum of GRF fluctuations (mini- 
mum determined by resolution) 
Total RMS amplitude of GRF fluctuations 
Relates ordered to isotropic random compo- 
nent 


Thermal electrons 






n c (r, <j>) = n c0 


Alternative constant density test model. 
By default, we use the NE2001 model of 
Cordcs & Lazioj |2002|), which peaks in the 
arms around 0.1cm -3 . 


Cosmic-ray electrons 


V 

h r 


3 

15 kpc 

0.25 
GcV s sr 


JcreOO = Jcre,® exp(-(r - R m )/h r ) 


Electron power spectrum power law index. See 
§0and §[33] 
Scale radius of CREs. 
See §[331 



Table 1. Table of modelling parameters as described in § [3] 



(i.e., no grid refinement with distance) and a simple 2D sim- 
ulation in the plane of the Galaxy only. We compute the 
observables at each of 512 pixels along the plane that cor- 
respond to the centres of the HEALPix pixels on the plane 
in a map at A/g^e = 128 map. We then smooth the result 
in the same way we smoothed the data as described in § 12.11 
and bin them into 64 pixels. 

4.2 MCMC with cosmomc 

Fo r the MCMC sampling , we use the tools already developed 
bv lLewis fc Bridle] (|2002r ) and publicly available as the pack- 
age COSMOMC. It consists of a Fortran sampling routine, into 
which one can insert a parametrized likelihood function, and 



an analysis tool, getdist, to read the resulting Markov chains 
and determine the parameters' mean and maximum likeli- 
hood values and correlation matrices. We use the COSMOMC 
sampler but have written our own analysis tool in order 
to add further functionality described below. The sampler 
has the ability to continuously update the proposal density 
based on the likelihood space so far sampled, which makes 
for a more efficient mapping out of the likelihood. 

This method is complicated in this case by the fact that 
we do not compute the model, i.e., the simulated emission 
and rotation measures, analytically. Rather, as described 
above, we simulate a Gaussian random field for the ran- 
dom component and apply the compression algorithm. The 
model we wish to compare to the data is then the expecta- 
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tion value of the emission profiles resulting from such simu- 
lations, and the corresponding variance among realisations 
determines the error. These are difficult to compute analyt- 
ically, so we instead use a brute-force approach of creating 
a set of N = 10 realisations and taking the mean. This 
must be done for each likelihood evaluation. (We have tried 
several different numbers of realisations at each sample and 
found ten to be the most efficient in terms of the trade-off 
between sample variance and run time.) 

For simplicity, we assume a Gaussian distribution of 
the data points, Xi, about the mean model value, fii, with a 
variance of Oi so that the likelihood is: 

- In (£) = In (v^a*) + (a* - /(2a*) (14) 

i 

(where the summation is over the pixels along the plane). 
The distributions are not exactly Gaussian, as can be seen 
in Fig. HJ] but the assumption appears to work reasonably 
well nonetheless. The differences come from the relation be- 
tween the Gaussian random field and the observables, which 
in the case of the synchrotron emission, is non-linear. The 
RMs of extragalactic sources start with an average value 
through the galaxy, dependent only on the coherent com- 
ponent of the magnetic field, that is then perturbed by 
a one-dimensional random walk due to the random mag- 
netic field component projected along the LOS (assuming 
the small-scale fluctuations are uncorrelated with the ther- 
mal electrons) . The random component then adds a variance 
proportional to the number of steps. (Because the random 
component has structure on many scale lengths, it is essen- 
tially a superposition of random walks with different sized 
steps.) For polarised intensity, the polarisation vector be- 
haves like a two-dimensional random walk, and the observ- 
able is only the length of the result. With only a GRF, the 
expectation value of either Stokes Q or U would be zero, but 
the expectation of P = \J Q 2 + U 2 is non-zero and adds to 
the coherent component. For total intensity, the emission of 
all components simply adds. 

Note that the variance at each longitude depends both 
on the position and on the model parameters and is deter- 
mined simply from the variance of the simulated realisations. 
There is an additional uncertainty in both the model mean, 
Hi, and its variance, <Ti, due to the limited number of reali- 
sations, so in the determination of x 2 , we actually use 

a 2 = (1 + 1/N)a1 (15) 

in the likelihood evaluation of Eq. 1141 

When analysing the chains, there is essentially an addi- 
tional "noise" term added to the likelihood space, since two 
evaluations at the same point in parameter space will re- 
turn slightly different likelihoods. This means that we need 
more samples to overcome the sample variance and to map 
accurately the shape of the likelihood space. It also means 
that we need to give a "temperature" parameter to the sam- 
pler. The latter is necessary to prevent the algorithm from 
falling into an unusually low \ 2 hole an d never getting out 
again. So rather than sampling from P oc exp— x 2 /2, we 
sample from P oc exp— x 2 /2T, where T is not a physical 
temperature but simply a way to adjust the sampler. This 
temperature is simply tuned to ensure that the acceptance 
rate of the sampler (i.e., the fraction of proposed samples 
that are accepted) remains roughly a third. (With a total 
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Figure 6. Distribution of data values compared to the mean. 
For 1000 simulations and 64 pixels along the plane, the histogram 
shows the difference between the data for a given simulation and 
at a given pixel, Xi, and the mean at that pixel over all simula- 
tions, rrii, normalised by the variance, cr;. Clearly, for total and 
polarised intensity, there is a tail on the positive side that biases 
the mean. The solid line shows a Gaussian of the equivalent mean 
and width, which is essentially what using a \ 2 assumes, while 
the dotted line shows a Gaussian that approximately fits the left 
half of the histogram. 



of 140 data points and corresponding \ 2 values, we find a 
temperature of 10 usually sufficient.) This temperature cor- 
rection is then reversed to analyse the results, correcting 
both the likelihoods as well as the weights (the number of 
steps in the chain where it remained at a given sample before 
it accepted the next). 

To simplify the analysis of the chains, we need to re- 
move this noise, and we do this simply by smoothing th e 
likelihood space (as, for example, done in lDick et al.| [2006'). 
For each sample, we replace the likelihood with the mean 
of the nearby likelihoods. We must also smooth the weights, 
which are by construction meant to be proportional to the 
likelihood. The result of the smoothing is a distribution of 
likelihood values approximating that seen in an essentially 
noiseless situation with otherwise identical parameters. We 
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smooth by binning the samples in likelihood space with a 
binsize of 1/15 times the range (after burn-in) for each pa- 
rameter. 

The number of samples required for the 
chains to converge (using COSMOMC's parameter 
MPI_Limit_Converge_Err = . 3) is between a few thousand 
for a 2D fit and a few hundred thousand for a 6D fit. 

In Sj5] we describe how we tested with simulated inputs 
to verify that the MCMC search returns the input param- 
eters and reasonable uncertainties. An example is shown in 
Fig.0 

4.3 Parameters and datasets 

One could, in theory, throw all data and all parameters into 
the MCMC analysis at once, but it would be surprising 
if such an approach yielded useful information. Degenera- 
cies among parameters make samplers highly inefficient, and 
when starting values are very far from the peak likelihood, 
the burn-in time to random-walk over to the peak can be 
very long. Instead, we make educated guesses (informed by 
trial and error) about what combinations of parameters and 
data will give good constraints and test them using sim- 
ulations. For the parameters we fix for the course of the 
analysis, we choose by hand values that roughly match the 
data, since these are not our primary interest. 

For example, as described above, the RMs are the best 
dataset to constrain the coherent magnetic field component. 
Particularly, the angle 0o describing the azimuthal orienta- 
tion of the spiral may have a strong effect on the position 
of the RM sign reversals, while the a n describe the relative 
magnitudes of the RM features. While these parameters do 
have an effect on the profiles of the synchrotron emission, 
given the galactic variance, this may not be detectable. So 
to determine these parameters, then, we use only the RM 
data. 

The synchrotron emission is particularly interesting for 
studying the different components of the random field. As 
described above, we can use the RMs to constrain the co- 
herent field and then use the synchrotron to constrain the 
parameters such as B TU1B and / or d, keeping cf>o and a n fixed 
to their best-fit values from the RM-only analysis. 

Table [1] lists all of the parameters that go into the mod- 
elling, and clearly studying all of them would be a huge 
undertaking. For this first paper, we focus on the most im- 
mediately interesting issues of the magnetic field reversals 
and the relative contribution of the coherent, ordered, and 
random components. 



5 TESTING 

5.1 Basic Parameters 

Figure [7]shows an example of how the MCMC method works 
with a simulated dataset. The results on the left use only 
RM data covering roughly half of the plane. Among the six 
parameters simultaneously fit for this test (all a n and (j>o), 
all return a mean, fj,i, that is within 2a of the input value, 
Ho, where a is the variance of the samples. 

On the right of Fig[7]is the result of the second step of 
fitting. We fix the amplitudes at the values found using the 



RM data only and then fit to all data the two parameters 
controlling the ratio of ordered to random magnetic fields in 
the arms. The returned values lie roughly within 1 a of the 
correct values. 

We have run a set of ten such simulations to verify that 
the estimates of / or d and Brms are correct and unbiased 
(i.e., | (m) -no\ < <t/vT0). 

5.2 Resolution 

Figure [8] shows a test of the importance of the resolution 
level in the simulations along the plane. Here, we test the im- 
plications of the simplifications we make for computational 
efficiency. 

The choice of the resolution of the GRF simulation is an 
important question. The GRF normalisation, e.g., its RMS 
over the whole galaxy box, is a function of the outer scale, 
the power law index, and the resolution. Though the smaller 
scales are increasingly irrelevant due to the power law's neg- 
ative spectral index, there is an effect of not including all rel- 
evant dynamical scales. We have verified that quadrupling 
the resolution appears to have only a small effect on the re- 
sults of the simple Galactic plane analysis, as seen in Fig. [8] 

More importantly, we simulate only a 2-dimensional 
GRF rather than a full 3-D galaxy, which speeds up the 
computation by almost a factor of ten. Though we are only 
looking along the Galactic plane, structure near the plane 
does enter the beam of each observation, particularly for vol- 
umes further away. To be fully correct, we would simulate 
a 3D galaxy and include contributions from high-resolution 
pixels just off the plane by using multiple HAMMURABI shells 
to simulate a finite instrument beam. The differences in the 
profiles are shown in Fig. [8] There is a visible difference in 
the total intensity profiles; the blue curve shows the correct 
profile, while the red is what we simulate. The effect of this 
simplification is to underestimate the RMS of the random 
component. If we perform the MCMC analysis to fit for the 
Bums and / or( j parameters, we find that the latter is unaf- 
fected but the former is underestimated by approximately 7 
per cent. We consider this acceptable given that our knowl- 
edge of the relative strengths of these components has not 
previously been anywhere near that accurate and that there 
are further limitations of our analysis discussed in § [6] 

5.3 Non-linear parameters 

Some of the parameters are fairly linearly related to the ob- 
servables. The arm amplitudes, for example, simply increase 
the RM in a given direction, and since the RMs are the best 
dataset to use to constrain them, we use only the RM data 
to determine the a n parameters (as seen in Fig. [7J| as well 
as the angle 0q. 

Parameters related to the synchrotron emission, partic- 
ularly those describing the random magnetic field compo- 
nent, are more clearly non-linear. Firstly, the emission itself 
depends on the field to the power of 7 as described in Eq. [1] 
Secondly, the observable polarised emission is affected by 
the random component much like a 2D random walk, whose 
mean distance from the starting point is then a function 
of the square root of the number of steps, N, and the step 
size. In the case of a power law GRF, we have essentially a 
superposition of random walks of many step sizes and Ns. 
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(c) a 2 (d) % (e) 

Figure 7. Example of MCMC results. The filled contours represent the mean likelihood showing the sampled region, while the white 
contour lines indicate the 68, 95, and 99.7 per cent confidence regions. The filled black diamonds give the correct input values. The 
examples on the left, (a) through (d), show the results for several of the six parameters describing the coherent magnetic field found 
fitting only to the RM data. On the right (e) are the results for the two parameters describing the amount of isotropic random and 
ordered components found using all the data. 



Furthermore, the variance used in the likelihood calcu- 
lation (Eq. I14[) depends on the particular location in param- 
eter space. As described above, this affects the "noise" due 
to the non-analytic way we compute the model. We have 
tested with a set of simulations to see how well the method 
works despite these shortcomings and found that it works 
surprisingly well. We fixed the other parameters and fit only 
to / or< j and -Brms, the two most important but perhaps most 
complicated and correlated parameters. Using both the total 
and polarised intensities does succeed in constraining these 
parameters. In all but one of the ten simulations, the result- 
ing best-fit positions were within 2-<r of the correct values 
and showing no bias. An example is shown in Fig. [3 



6 RESULTS 

Our principle aim is to attempt to constrain the ratios of 
the three components of the Galactic magnetic field: the co- 
herent component, the isotropic random component, and the 
ordered component. With the three complementary datasets 
of total synchrotron emission, polarised synchrotron emis- 
sion, and Faraday rotation measures, we have demonstrated 
that these three components can be disentangled using rea- 
sonable models of the magnetised ISM. 

6.1 RM fits 

We begin by constraining the large scale, coherent field. In 
other galaxies, the magnetic field usually shows spiral struc- 
ture, sometimes azimuthally symmetric, sometimes not. We 
should not assume that the magnetic spiral follows the DIG 



spiral; the field strength may or may not be higher in DIG 
arms, and even the spiral pitch angles may not be the same. 
But it is not our primary aim in this work to determine the 
global morphology of the coherent magnetic field. We use the 
RM data largely to constrain the amplitude of the emission 
due to the coherent field. Though its morphology may not 
follow precisely the spiral model we have fit to the data, it 
is enough to limit the amount of synchrotron emission from 
that component, allowing us to then examine that from the 
random and ordered components independently. 

Using the methodology described in § [4] we fit a six pa- 
rameter spiral model to the RM data alone to constrain the 
values of the following parameters: the orientation angle, 4>o, 
the amplitudes ao to 0,3 for each of the four arms, and 04 
for the "molecular ring". The amplitudes a„ then give the 
relative strengths of each feature, which combine with Co to 
give the maximum field strength along the ridge. The rela- 
tive amplitudes, particularly their signs that determine the 
direction of the field in each arm (clockwise or counterclock- 
wise when viewed from the north Galactic pole) determine 
the features seen in the RM data as in Figs. [5] or [S] 

The parameter fit results are shown in Table [2] This 
part of the analysis does not give much that is n ew compared 
to what has been done bv lBrown et al.l |2007l ) and others. 
We do clearly see the magnetic field reversals, visible in the 
RM data, in the general direction of Scutum-Crux and near 
the tangent to the proposed "molecular ring" at negative 
longitudes. This is visible in comparing Figures [2] and [4] We 
also confirm the field strength peaking at ~ 2 /iG in the arm 
ridges. It is interesting to note that, though the field is free 
to fit the azimuthal rotation </>o, the result lies very close 
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Figure 8. Testing the effects of resolution and dimension of GRF 
simulation as described in § 15,21 The solid lines are the mean mod- 
els while the dotted lines indicate the RMS. The difference be- 
tween our 2D analysis and a better but more time-consuming 3D 
method is apparent only in the strength of the step-features cor- 
responding to spiral arm tangents where the small-scale random 
component is strongest. These differences are within the uncer- 
tainties but are systematic and introduce a small bias discussed 
in §E3 

to the DIG spiral arms in the NE2001 model, whether that 
model is assumed for the thermal electrons or not. 

The value of ao is not constrained, because as seen in 
Fig. [4l no sightlines are really tangent to the arm. It con- 
tributes relatively little in the region where we have RM 
data. 

6.2 Synchrotron fits 

Our main aim, however, is to attempt to quantify the rel- 
ative amounts of coherent, random, and ordered magnetic 
field components. As found by Broadbent at al., the mag- 
netic field must be amplified in the arms. Since the cos- 
mic rays are likely smooth on these scales, the step-like fea- 
tures in the synchrotron profile are probably due to the en- 
hancement of the perpendicular component of the magnetic 
field in the arms. Though the coherent component could be 
scaled to produce such structures, it would also require an 
extremely low (n c w 10~ 3 cm" 3 ) electron density in order 
to be consistent with the RM data. 

The fit results in Table [5] for the parameters Arms and 
/ or d are simply due to the relative strengths of the total 
and polarised emission along the plane. The resulting best- 



fit profiles and residuals are shown in Fig. [TO] These profiles 
are not a perfect match to the total intensity data, but they 
agree remarkably well in most places. The most significant 
deviation is where there is a strong thermal feature around 
longitude -70° to -80°, which is the region that dominates 
the high % 2 ■ There is also a mismatch in the profile around 
longitude 60° which bears further investigating but appears 
to be related to the Sagittarius arm. Recall from § 12. II that 
the innermost region is not included in the fit. Further work 
will be necessary to determine the magnetic field structure in 
the Galactic centre region; the model we use, particularly the 
central part of the profile shown in Fig. [3] is ill-constrained. 

Figure [10] shows that the polarised synchrotron profiles 
are only fit to first order. The x 2 values are not bad due to 
the galactic variance, but the residuals show clear system- 
atic differences where the profile shape is wrong. More work 
will be needed to determine what this means for the model, 
particularly the spatial variation of the magnetic field com- 
ponents. 

6.3 Component ratios 

From the results of our fitting the parameters cj>o, a n , B rms , 
and /ord, we learn the following: 

(i) The rotation measures allow us to constrain the co- 
herent magnetic field parameters a n . As shown in Fig. |H 
they alone cannot distinguish between the NE2001 spiral 
model for the thermal electron density and a simple uni- 
form model. Not shown on that figure are the synchrotron 
profiles, which are effectively indistinguishable. The reason 
is that the synchrotron profile is apparently dominated by 
the random and ordered components. The RMs, regardless 
of the thermal electron model, have constrained the coher- 
ent field to an amplitude that is fairly low compared to the 
synchrotron profile, assuming that the cosmic ray electron 
density is roughly correct. 

Therefore, our results suggest that the coherent field in 
the arms peaks around Boa n (Co + 1) = 1 — 3 fiG in the 
arms, varying from arm to arm with an average of 2 /j,G. 

(ii) The isotropic random field component peaks around 

(A 2 so) 1/2 = (Co + 1)S RMS = 4.2 /iG (see Eq.[U). 

(iii) The ordered field component peaks slightly lower at 
3.3 mG. 

(iv) The fractions of the energy density in coherent, ran- 
dom, and ordered field components (oc (B 2 )) are then 
roughly 1:5:3, respectively. The uncertainties are discussed 
further in § 16.61 Note that this model implies a total local 
field strength of roughly 3 [iG, and a mean field strength in 
the inner 10 kpc of the Galaxy of 3.5 /xG. 



6.4 Arm strength 

In the basic analysis, the arm strengths, a n , determine the 
coherent field strength in each arm (as well as direction) 
but do not affect the random or ordered components. For 
simplicity, these components are simply scaled by p c as de- 
scribed in Table [JJ which is only a function of the distance 
from the ridge. We can also, however, adjust the strength of 
these components by the same factors by multiplying Co by 
\a n \. This makes the contrast in each arm proportional to 
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00 


70 ±6.2 


70 ± 7.8 


do 


1.65 ± 0.38 


-0.15 ±0.24 


ai 


0.61 ± 0.18 


0.33 ± 0.15 


a 2 


-1.04 ±0.21 


-0.87 ±0.20 


A3 


1.26 ± 0.64 


0.97 ±0.46 


4 


-1.00 ±0.35 


-1.33 ±0.43 


Fits to all data with above NE2001 


n e values fixed 




constant compression 


scaled 


-Brms 


2.1 ± 0.03 


2.1 ±0.04 


/ord 


1.9 ±0.14 


1.5 ± 0.16 


xi/N 

XlJN 
X 2 KM /N 


3.7 
0.8 
1.8 


3.5 
1.0 
1.6 



Table 2. Top: results of fitting RMs with the simple spiral 
model parameters and comparing constant thermal electron den- 
sity model to NE2001. Bottom: results of fixing the 4>o and a n 
parameters to the best fit values show on top using NE2001 and 
all others to the "defaults" listed in Table [T] and then fitting 
only the / or[ j and Brms parameters to all datasets. As described 
in § 16.41 the "constant compression" have the same amplitude 
of random components in each arm while the "scaled" have the 
factors of a n applied. 



the coherent field strength of that arm rather than constant. 
(It also implies that the ratio of the ordered to isotropic ran- 
dom fields changes from arm to arm.) 

A comparison of the best-fit profiles is shown in Fig. 1101 
The physical motivation for scaling each arm strength in- 
dependently is unclear, but the result including the extra 
scaling appears to fit the synchrotron profile slightly better 
in some places and slightly worse in others. 



6.5 Arm/inter-arm Contrast 

In the basic analysis presented above, we have fixed the value 
of the parameter Co that defines the contrast between the 
arm ridges and inter-arm regions. We can instead allow this 
parameter to vary as well. The combination of Co and Brms 
controls the relative strength of the arm versus the inter-arm 
region for the random component (as Co and Bo or a n does 
for the coherent component). 

These three parameters, however, are somewhat degen- 
erate as shown in Fig. 1111 The data do not easily distinguish 
between a stronger Brms with a weaker contrast (Co) or vice 
versa. This is an example of one of many degeneracies in the 
parameters listed in Table [T] But our basic result comparing 
the strengths of the components in the arms is not signifi- 
cantly affected. 



6.6 Discussion 

We have used the three complementary datasets of total syn- 
chrotron emission, polarised synchrotron emission, and ro- 
tation measure to study the three components of the Galac- 
tic magnetic field, namely the coherent, random, and or- 
dered fields. We find that we can fit the profiles of these 
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Figure 9. Top, (a) through (d): selected results of fitting the 
coherent field parameters to the RM data only using the NE2001 
thermal electron density model. Bottom, (e): results of fixing the 
coherent field parameters to the values found in the RM analysis, 
and then fitting the two parameters controlling the small-scale, 
irregular components to all the data. The means and uncertainties 
are given in Table [2] Compare to Fig. [7] 



datasets along the plane with a model where the compo- 
nents peak in magnetic arm ridges and where the coherent 
component contributes roughly 10 per cent to the energy 
density (<^B 2 )), while the random and ordered components 
contribute roughly 50 and 40 per cent, respectively. 

This is a first attempt to constrain models with many 
parameters and degeneracies. If the estimates for the CRE 
and thermal electron densities are even approximately right, 
our method gives a fairly robust constraint on the relative 
strengths of the three field components. We should, how- 
ever, quantify the impact of such assumptions and their un- 
certainties: 

• There is a degeneracy between the thermal electron 
density in the arms and the amplitude of the coherent mag- 
netic field. A factor of two uncertainty in the mean electron 
density in the arms implies a factor of two uncertainty in 
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Figure 11. MCMC results showing the degeneracies described in £j 16.51 among the three parameters: Bp,MSi /ordi an d Co- 



the coherent field strength. If the electron density is a factor 
of two smaller than our model for a given arm, the coherent 
field is then a factor of two larger to reproduce the same RM 
profile. This then requires a drop in B r ms to reproduce the 
same synchrotron total intensity profiles. The ratios of en- 
ergy densities of the components (coherent :random:ordered) 
would go from roughly 1:5:3 to 6:3:1. 

• Similarly, in using the NE2001 model for a smoothly 
distributed, average electron density, we have assumed a 
relatively homogeneous DIG, while clearly the density of 
free electrons varies many orders of magnitude from HII re - 
gions down to molecular clouds. Bcrkhuiis en et al.l (2006) 
and others have studied it's dumpiness as described by 
a filling factor. The LOS filling factor is often defined as 
/ = (rig) /n 2 ., where the brackets denote the average along 
the entire LOS while n c is the average value in the individ- 
ual clumps; it then also represents the fraction of the LOS 
that falls within the clumps. This can also be defined as 
/ = (n e ) 2 / {nl) = DM 2 /{EM D) (where D is the total 
LOS distance, EM and DM are the emission and disper- 
sion measures, respectively). We note that the factor de- 
rived from comparing EM and DM oversimplifies the fact 
that EM and DM essentially trace different physical compo- 
nents of the ISM. 

The NE2001 is largely based on DM, which is linear in 
n e and path length, and so any factor of / 7^ 1 cancels. If 
the fluctuations in n e are uncorrelated with the magnetic 
field, then RM is similarly unaffected on average. It is un- 
clear how correlated the field is with the electron density, 
however, so the impact of this approximation is difficult to 
judge. If the field fluctuations along the LOS are equally 
likely to be toward or away from the observer in high den- 
sity clumps, then even if the strength of the fluctuations is 
correlated with n e , this will not affect the predicted mean 
RM profile, only its variance. If, however, the correlated field 
fluctuations retain a coherent direction, then the observed 
RM would be larger than what we have modelled. Properly 
taking this into account would then require a weaker field 
to reproduce the same RM data. In future, we may be able 
to plug in the results of magneto-hydrodynamic (MHD) tur- 
bulence simulations into the tools we have described to test 
this explicitly. 

• There is a degeneracy between the magnetic field 
strength and the cosmic ray electron density. We assume 
a value of Jcre,® = 0.25 (GeVm 2 ssr) at the position 
of the Sun. If the uncertainty in this number is a factor of 



two (which it may well be; see, e.g.. lStrong et af1l2004l ) and 
assuming B co h is constrained by the RM data, then we can 
solve for the other components. The ratios of energy den- 
sities would then vary from 0.6:5:4 (CREs halved) to 3:5:2 
(CREs doubled). 

The recent re sults from the Fermi 7-ray telescope 
l|Abdo et al.ll2009D seem to indic ate a density that is some- 
what lower than what we use from lStrong et all (|2007l ) based 
on EGRET data. A lower CRE density implies a larger ran- 
dom magnetic field component to re produce the obse rved 
synchrotron emission. (The work of 1 un et al.l (|2008l ) use 
an even higher density of Jcre,® = 0.4 (GeVm 2 ssr) , 
which partly explains how they can reproduce roughly the 
observed amount of emission with only an isotropic random 
component and no ordered component.) 

• We have assumed a spatially smooth cosmic-ray elec- 
tron distribution with a simple power-law power spectrum 
that is the same throughout the Galaxy and has index p — 3. 
A spatial variation in the spectrum or deviations of the spec- 
trum from the assumed power law would also affect our anal- 
ysis. 

As an example, we can imagine that the CRE spatial and 
spectral distribution varied between spiral arm and inter- 
arm regions due to the higher density of CRE accelera- 
tion sites (e.g., supernova r emnants) in the arms (see, e.g., 
ICase fc Bhattacharval 0*9961 ). This would have a similar ef- 
fect as the global uncertainty discussed above. 

Furthermore, in assuming this power law distribution of 
electrons, we assume that the spectrum of synchrotron emis- 
sion in brightness temperature also follows a power law with 
/3 = — ( p + 3)/2 = —3 . Anal yses such as iPlatania et ail 
( 1998) or iGiardino et~ai1 (|2002T l show that this index in fact 
varies across the sky. Determining this index on the plane 
over a broad frequency range is complicated by the variety of 
emission mechanisms that contribute to total intensity and 
by the Faraday effects that contribute to polarised intensity 
at low frequencies. If our assumption of ft — — 3 is too steep, 
our model will be underpredicting the polarised emission at 
23 GHz from the coherent magnetic field component, for 

example. This is difficult to quantify, however, since both 

p+i 

the synchrotron emissivity (as B~) and the polarisation 



fraction (II — 



g+i 



depend on p. 



p+7/3' 

• There are uncertainties in the zero-levels of the syn- 
chrotron datasets we used. These uncertainties are relatively 
small (e.g., a few K for the 408 MHz map), usually much 
less than the galactic variance, and will be most significant 
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Figure 10. Best fit models from MCMC analysis as described 
in § 16,11 and 16,21 The three observables are shown from the top 
as synchrotron total intensity at 408 MHz, synchrotron polarised 
intensity at 23 GHz, and rotation measure. The solid blue and 
magenta lines are the mean models while the dotted lines indicate 
the RMS, or galactic variance, and in thick green are the data. 
For each observable, beneath the plotted profile are the residuals, 
i.e., (data — model) / 'error , compared to dashed lines at ±1 and 
±2 a. In magenta is the result where all arms have equal random 
and ordered compression, while in blue is the result where those 
components are also scaled by (lit (IS described in § 16.41 



where the synchrotron emission profile is lowest toward the 
Galactic anticentre. For our analysis of the magnetic field 
component ratios in the spiral arm peaks, this uncertainty 
is certainly not significant. It will, however, affect the pa- 
rameters describing the strength of the field and the density 
of CREs in the outer region of the galaxy. As we have dis- 
cussed above, these two distributions are degenerate and 
not constrained by our analysis. In future work, when we 
add dust emission to the analysis to break the degeneracy, 
we will also have to address the zero-levels more accurately. 
For this work, however, we consider the impact of the un- 
certainty minimal. 

• We have assumed that the large-scale coherent field and 



the small-scale field components are all strongest in the same 
regions, the magnetic spiral arm ridges. There is some ev- 
idence in one galaxy that the coherent field, in contrast to 
the random field, may be stronger in the inter-arm regions 
|Beck fc Hoerneal 19961 ). Further work will be needed to de- 
termine whether we can distinguish these possibilities. 

Our results are the ratios of the field components' en- 
ergy densities that we measure as roughly 1:5:4 (coher- 
ent:random:ordered). The above uncertainties most signif- 
icantly affect the fraction of magnetic field energy in the 
coherent component. The rough estimates above then imply 
that the coherent component could make up as little as 5 
percent of the total energy in the field, or as much as 30 
per cent given an uncertainty of a factor of 2 in the CRE 
density, or as much as 60 per cent if the thermal electrons 
are half the assumed value. The latter is certainly possible 
for a small region, but on the whole, pulsar distance mea- 
sures averaged over the galaxy are unlikely to allow such 
a large discrepancy. As mentioned above, Fermi estimates 
of the CRE density are even lower than what we assumed, 
implying even larger random and ordered components. 

There have been other efforts to estimate the amount 
of turbulence in the magnetic field. But one must be care- 
ful to understand what it is that is being measured, as 
these works generally distinguish only between a "regu- 
lar" field and a "turbulent" field. Depending on the observ- 
able, what we define as the o rder ed field may count as e i- 
ther. iRand fc Kulkarnil (|l989T ) and lOhno fc Shibatal l|l993h . 
for example, both look at RM data and measure a tur- 
bulent magnetic field component of 4 — 6 pG that corre- 
sponds to t he sum of our or d ered and isotropic random 
components. ISchnitzeler et af] l|2007h give an upper limit 
of -Bturb/B rcgi x ~ 2 based on synchrotron emission. Their 
"regular" component is, when only looking at synchrotron 
emissi on, the total of our c oherent and ordered compo- 
nents. [Haverkgrn^^Tjil] {2004), by contrast, find a random 
component smaller than the "regular" component, though 
this is in two relatively sm all fields that may be probing 
smaller-scale structure only. iMiville-Deschenes et alj (|20Qg ) 
find Bturb/Brcg,± = 0.57 using synchrotron polarisation in 
WMAP data, so again, the "regular" component includes 
what we define as ordered. Each of these studies uses differ- 
ent observables, different regions of the Galaxy, and different 
assumptions. It is thus difficult to compare their results with 
each other or with our estimate, but with the exception of 
the Haverkorn et al. result, they are roughly consistent. 

What we can conclude, however, is that the coherent 
field is likely a relatively small fraction of the total and that 
a model including only an isotropic random component and 
a coherent field is incomplete. To reproduce all of the ob- 
servables, particularly the polarised emission, requires an 
ordered field component as well. 



7 CONCLUSIONS 

We have outlined a method to simulate observables such as 
synchrotron emission and Faraday rotation measure based 
on parametrised models of the physical components of the 
Galaxy and to use these to study the complicated parameter 
space in an MCMC analysis. In particular, we describe how 
to deal with the galactic variance due to the fact that one of 
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the principle physical components is stochastic. These tools 
can be used to model the global properties of the Galaxy's 
magnetic field in 3D or to study the turbulent component 
of the magneto-ionic medium in detail without the need for 
as many simplifying assumptions. 

To demonstrate the method's utility, we have used the 
three complementary datasets of total synchrotron emission, 
polarised synchrotron emission, and rotation measure to try 
to constrain the relative contributions of the coherent, ran- 
dom, and ordered components of the magnetic field. We 
find the relative energy densities of these components to 
be roughly 1:5:3, respectively, at their peaks in a magnetic 
spiral arm model. If our assumed CRE and thermal elec- 
tron distributions are roughly correct, then this implies arm 
field strengths of 2, 4, and 3 fj,G, respectively. This partic- 
ular analysis is limited by several simplifying assumptions, 
but we consider the result a first step and proof-of-concept 
showing how to tackle this complicated problem. 

In future work, we will address some of the limita- 
tions and degeneracies discussed above. More importantly, 
we look forward to the prospect of additional data. In the 
next few years, we will have the C-Band All Sky Survey 
(C-BASSj3 at 5 GHz giving a much more sensitive mea- 
surement of the polarised synchrotron emission in the radio 
bands. These data will help us to study the spectral variation 
of the synchrotron emission, and in turn the distribution of 
cosmic ray electrons. 

The Planck satellite (jThe Planck Collaboration! l200d ) 
will also provide a higher sensitivity polarised synchrotron 
sky map in the microwave bands when combined with con- 
tinuing observations by WMAP. Its many bands will also 
contribute, along with C-BASS and WMAP, to a much more 
accurate separation of the thermal emission than is currently 
possible. This will improve our fitting of the step features in 
the synchrotron emission along the plane where the thermal 
emission is a problem. Furthermore, the High Frequency In- 
strument on Planck will give us a much better map of the 
polarised dust emission, a completely independent tracer of 
the magnetic fields that we have not used in this work. 

The GALFACTS Q survey at Arecibo Observatory will 
give us the vital coverage of the northern side of the Galac- 
tic centre in rotation measures (see Fig. |2j. This will be 
crucial in distinguishing the various models of the magnetic 
field reversals currently so controversial. And that in turn 
will inform theories such as dynamo models for the origin, 
amplification, and evolution of the coherent field. 

The Australian Square Kilometre Array Pathfinder 
(ASKAP) 

will significantly improve our southern sky 
coverage of rotation measures. Even better will be the 
Squ are Kilometre Arra^Q itself. Ijohnston fc et al l (|2008l ) 
and iGaensler et alj (|2004 ) describe how these will signifi- 
cantly advance magnetic field studies. Not only will they in- 
form our knowledge of our own galaxy but allow us to make 
similar observations of external galaxies using background 
polarised sources. 



4 http : / /www . astro . caltech . edu/cbass/ 

5 http : //www . ucalgary . ca/r as/GALFACTS/ 

6 http: / /www. atnf . csiro . au/projects/askap/ 

7 http://www.skatelescope.org/ 
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